On the detectability of non-trivial topologies 
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We explore the main physical processes which potentially affect the topological signal in the 
Cosmic Microwave Background (CMB) for a range of toroidal universes. We consider specifically 
reionisation, the integrated Sachs- Wolfe (ISW) effect, the size of the causal horizon, topological 
defects and primordial gravitational waves. We use three estimators: the information content, the 
S/N statistic and the Bayesian evidence. While reionisation has nearly no effect on the estimators, 
we show that taking into account the ISW strongly decreases our ability to detect the topological 
signal. We also study the impact of varying the relevant cosmological parameters within the 2a 
ranges allowed by present data. We find that only Qa, which influences both ISW and the size of 
the causal horizon, significantly alters the detection for all three estimators considered here. 

PACS numbers: 
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I. INTRODUCTION 

One of the fundamental questions of cosmology is 
whether the universe is infinite or finite, a question which 
is especially relevant in the context of string theory where 
several dimensions may be compact. With the move 
towards precision cosmology, it has become possible to 
study the topology of the universe up to the causal hori- 
zon which fundamentally limits our view of the world. 

Many tests have recently been developed to detect and 
constrain the topology of the universe, see Refs. [l|,[2j and 
references therein. Observationally a non-trivial topol- 
ogy can be detected directly with the observed map of 
CMB temperature and polarisation fluctuations through 
matched and correlated circles in the CMB d, 0, H, IE 
0, The topology can also be searched for in the 
expansion of the CMB map in terms of spherical har- 
monics, T(x) — ai m Y£ m (x), where x are the pixels 
and both the temperature field and the ag m are random 
variables. The presence of a non-trivial topology intro- 
duces preferred directions, which breaks global isotropy 
and induces correlations between off-diagonal elements 
of the two-point correlation matrix (a£ m a|, m ,). Methods 
to find such off-diagonal correlations and to put optimal 
constraints on the size of the universe were proposed in 
Ref. 0. They are based on matching the measured corre- 
lations to a given correlation matrix through a correlation 
coefficient, as well as a maximum likelihood method and 
Bayesian model comparison. 

So far the focus was purely on trying to detect any 
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non-trivial topologies. However, to determine reliable 
limits, it is also necessary to understand all effects which 
can affect the detection. For this purpose, we study the 
impact of changing the cosmological parameters in the 
ranges allowed by present data. This should be seen as 
a step towards a rigorous search for the presence of signs 
of non-trivial topologies in the CMB. Eventually, such 
an analysis should account for other contributions like 
foreground emissions as well as sky cuts. 

One additional reason to perform such a detailed anal- 
ysis comes from statistics. The potentially most direct 
way to quantify the probability that our universe is infi- 
nite is to compute the model probability with Bayesian 
evidence. To do this, one needs to integrate over all pa- 
rameters, which means that one needs to know the be- 
haviour of the likelihood for a topology under changes of 
the parameters. 

This paper is organised as follows: we start by re- 
viewing our notation and the likelihood-based probes of 
topology. We then study the main physical effects that 
are relevant for topology studies. Finally we discuss how 
changes in the cosmological parameters affect our detec- 
tion ability, and what this means for realistic universes 
where the parameters are not known with absolute cer- 
tainty. 



II. METHODOLOGY 

In this section, all the simulations (as described in 
Refs. [H EU) are based on a flat ACDM model with 
Qa = 0.7, a reduced Hubble parameter of ft. = O.OIifo x 
(1 km s- 1 Mpc" 1 )- 1 = 0.67, a Harr ison-Zcl'dovich ini- 
tial power spectrum (ng = I) and a baryon density of 
ilbh 2 = 0.019. With this choice of cosmological param- 
eters we find a Hubble radius c/Hq w 4.8 Gpc, while 
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the radius of the particle horizon is i?h ~ 15.6 Gpc. In 
this paper we limit ourselves to cubic tori, denoted as 
T[X,X,X] where X is the size of the fundamental domain, 
in units of the Hubble radius. As an example, T[4,4,4] 
is a cubic torus of size (19.3 Gpc) 3 . The volume of such 
a torus is nearly half that of the observable universe. 
The diameter of the particle horizon is about 6.5c/Hq 
(cf. section UlI A[) . We will focus on three different sizes, 
T[2,2,2], T[4,4,4] and T[6,6,6]. We combine the £ and m 
indices of the spherical harmonic coefficients to a single 
index s = £(£+1) + m and mix both notations frequently. 

In this paper we use the likelihood based methods of 
Ref. @, which we review here quickly. The CMB likeli- 
hood is taken to be a multivariate Gaussian distribution 
with covariance matrix B ss > — (a*a s i) and zero mean, 
(a s ) = 0. This covariance matrix provides the full statis- 
tical description of the CMB for a given topology. The 
likelihood function then is 



To improve the stability of the algorithm against mis- 
estimates of the power spectrum we whiten the ag m by 
dividing out the power spectrum 
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The power spectrum of the ae m has to be estimated from 
the data, which changes their distribution. We found in 
Ref. @ that this tends to increase the sensitivity of the 
estimators discussed below. 



C(B)=p(a tm \B)(x 



1 



exp 



4e. 



IB-las 



(1) 

where \B\ is the determinant of the matrix B. Any further 
model assumptions like the values of the cosmological 
parameters arc implicitly included in the choice of B. 

Often it is easier to consider the logarithm of the like- 
lihood, ln£ = -l/2(ln|B| + 

X 2 + s max ln(27r)) where we 

have defined 



X 2 = £<S ss !a s 



(2) 



Writing (a*a s >) — A ss > for the expectation value of the 
two point function of the "sky" at m , we can compute the 
expectation value of the x 2 > 



(3) 



where we have used the hermeticity of the correlation 
matrices. For the two special cases of an infinite universe 
and one for which A = B (ie. the measured ae m are 
drawn from the template covariance matrix B) , we find 



<X 2 )oo = 

S 

(x 2 )b = tr(l) = s m8 
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As the ai m are Gaussian random variables, we expect 
to find that \ 2 is distributed with a x 2 -like distribution. 
The general expression for the variance is rather cumber- 
some, but for the two special cases we find 
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A. Information content 

As discussed in Ref. Q it is possible to study directly 
the information content in the CMB sky. For Gaussian 
fluctuations, the infinite universe contains a maximum of 
information (maximal entropy), while a finite ("small") 
universe is less random. Mathematically this is expressed 
through the non-zero correlations in the off-diagonal part 
of the two-point function which make it easier to predict 
the temperature distribution (at least in theory). 

A commonly used measure of information content is 
given by the Kullback-Leibler (KL) divergence between 
two probability distributions p and q, 



p\og [ - 
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This quantity measures the relative information entropy 
between the two distributions. It is often used in infor- 
mation theory. For example, if information distributed 
according to q is transmitted with a coding based on p 
then DKh{p\\q) describes how much bandwidth is wasted 
because of the wrong coding (the base of the logarithm 
is equivalent to a choice of units in which the KL di- 
vergence is expressed). It is also possible to view the 
KL divergence as a regularised entropy for continuous 
random variables, as it satisfies the Gibbs inequality, 
DxL{p\\q) > with equality only for p = q. Here we 
use the latter interpretation. 

In our case we are dealing with Gaussian probability 
distributions characterised by zero mean and covariance 
matrices A and B, for which the KL divergence can be 
evaluated directly, 



log |^| -tr(l-B-U)] 
(12) 



Unfortunately the KL divergence is not symmetric in its 
arguments, making its interpretation somewhat difficult. 
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It is also not rotationally invariant, so that it depends 
on the relative orientation of the coordinate systems in 
which the correlation matrices are expressed. We will 
therefore only consider a special case where we take p{B) 
to be the distribution of an infinite universe. In this case 
the KL divergence is the relative entropy to the topology 
described by the correlation matrix A. This is both ro- 
tationally invariant (as the infinite universe is isotropic) 
and has a well-defined interpretation. It is also very easy 
to compute: For the infinite universe B = B^ 1 = 1. The 
whitening enforces tr(„4) — s max = tr(f), so that the 
relative entropy of B is just 



H(A) = D KL (A\\1) 



1 



log |-4| . 



(13) 



This quantity describes how much less information is 
present in CMB data of the (easier to predict) finite uni- 
verse, as compared to the information of the CMB of an 
infinite universe. This gives an upper limit on the chance 
to detect a given topology. 

We plot in Fig. [1] both information content H as a 
function of £ max (left panel) as well as the increase of H 
per I (right panel) for the 3 topological models T[2,2,2] 
(red solid line), T[4,4,4] (blue dashed line) and T[6,6,6] 
(green dotted line) respectively. The information content 
decreases as expected for topologies closer to an infinite 
universe (T [6,6,6] in our case). 
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FIG. 1: The information content as a function of the size of 
the fundamental domain and of £ m a X : The left panel shows 
H and right panel the increase per £, AH, for (from the top 
downwards) T[2,2,2] (red solid line), T[4,4,4] (blue dashed 
line) and T[6,6,6] (green dotted line). 



For the first analysis we compute the correlation be- 
tween an observed map and a template and we search the 
global minimum of \ 2 over all rotations of the ai m (corre- 
sponding to relative orientations of the observed map and 
template). In order to quantify possible false detections 
of the topological signal, we give an estimate of the null 
hypothesis, that is the possibility for an infinite universe 
to be taken for a universe with a non trivial topology. 
We use the following procedure: For each given topol- 
ogy (T[2,2,2], T[4,4,4], T[6,6,6]) we fix the template and 
generate several thousand random maps both from the 
chosen template and for an infinite universe. We then use 
the mean and standard deviation for both distributions 
to compute a signal to noise statistic, 



S/N(B, X 2 ) = 



\(x 2 )qq-(x 2 )b\ 



<X 2 )\ 



(14) 



This corresponds to the distance between the means in 
units of the errors added in quadrature. We show results 
in Fig. [2l with the same colours/line styles as in Fig. Q] 
Again, we see that larger universes are more difficult to 
detect. The T[6,6,6] universe is barely detectable, while 
the others can be detect easily for £ max > 16. 
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FIG. 2: S/N for (from the top downwards) T[2,2,2] (red solid 
line), T[4,4,4] (blue dashed line) and T[6,6,6] (green dotted 
line) . 



B. Measuring topology 

While the previous subsection deals with a purely the- 
oretical measure of detectability, namely the information 
content in a given map or matrix, we apply in this sec- 
tion two actual detection algorithms discussed in Ref. @ . 
For both we use likelihood, Eq. (fTJ), or its logarithm % 2 , 
Eq. {5}. We exploit the fact that the the likelihood will 
be higher if observed map and template agree. 



The above test is somewhat ad-hoc, so that we com- 
pute additionally the probability that the ae m are actu- 
ally drawn from the distribution described by the correla- 
tion matrix B. This is often called the Bayesian evidence 
in astrophysical applications. Formally, the probability 
of a model (specified here via B) is p(B\ae m ), which we 
can connect to the likelihood, Eq. (fl]) through Bayes the- 
orem, 

P (B\a em ) = p{a lm \B)^-r cx C(B)p(B). (15) 
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We will be interested in ratios of the these probabilities 
for different B (describing different topologies). As we 
have no a priori knowledge of the topology, we will take 
p(B) to be constant. However, the numerical value of 
C(B) depends on the orientation of the template, so that 
we need to marginalise over the angles (see Ref. Q for 
more details, and [HI, 0, E3 for other examples of the 
use of Bayesian evidence in topology). Unfortunately 
the likelihood is very strongly peaked in the vicinity of 
correct alignments. The numerical integration is difficult, 
even using adaptive schemes. 

Having computed p(B\ae m ) we define a relative evi- 
dence for a topology as the ratio of this probability to 
the one for an infinite universe, 



Alog_E(S) = logp(B\a tm ) - logp(l\ae m ). 



(16) 



This gives the relative probability that the ai m are due 
to a universe described by the topology and cosmology 
encoded in B, compared to the probability that they are 
due to an infinite universe. Fig. [3] shows the strong de- 
pendence of the evidence on the size of fundamental do- 
main, the figure looks somewhat similar to the one of H. 
One can further eliminate the dependence on the cos- 
mological parameters by marginalising (integrating) over 
them as well. One of the motivations for this paper is 
to study the dependence of A log i? (23) on the cosmolog- 
ical parameters and to understand which integrals are 
necessary. Finally, by also marginalising over a class of 
topologies (e.g. the X of cubic tori T[X,X,X]) and using 
the actually measured ai m we can compute the probabil- 
ity that our universe has that kind of topology. 
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FIG. 3: Evidence for (from the top downwards) T[2,2,2] (red 
solid line), T[4,4,4] (blue dashed line) and T[6,6,6] (green dot- 
ted line) , compared to the atm of an infinite universe. 



III. PHYSICAL EFFECTS AT LARGE SCALES 

The shape of the universe leaves specific imprints on 
large scale CMB anisotropics at multipoles lower than 
£ ~ 60 that can be detected via the correlation matrix 
of the ai m coefficients. Such scales suffer many prob- 
lems. Firstly, it is at these scales that cosmic variance 
ACg = y/2/ '{21+ l)Ci is the largest. Then, at those 
scales several physical processes affect the CMB signal 
which may change our ability to detect non-trivial topolo- 
gies and to measure their properties. In the following we 
will study the main effects, namely the reionisation and 
the integrated Sachs- Wolfe (ISW) effect, in addition to 
the relative size of the particle horizon as well as tensor 
modes and defects. In the next section, we will vary the 
relevant cosmological parameters and study their influ- 
ence on the detection of the topological signal. 



A. Relative size of fundamental domain and 
particle horizon 

As discussed in the introduction, we fix the size of fun- 
damental domain in units of c/Hq. Another key scale is 
the distance to the last scattering surface which is very 
close to the size of the particle horizon. In the following, 
we will not discriminate between these two. As light- 
rays move on geodesies with ds 2 = we can compute the 
particle horizon today in comoving coordinates from 
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The quantity relevant to us is the diameter of the par- 
ticle horizon D (the "size of the causal domain"), which, 
in units of c/Hq is 



D = 



2a Q H Q r] Q 



= 2 



dx (u> r h 2 + Q m x + Q^x 4 ^ 



-1/2 



where the new integration variable is the normalized scale 
factor x = <x/ao, and we have assumed flat spacelike sec- 
tions and dark energy in the form of a pure cosmological 
constant. The energy density in radiation u> r = £l r h 2 
is fixed very precisely by the temperature of the CMB 
(assuming massless neutrini) to 4.2 x 10~ 5 . Its pres- 
ence changes the value of the above integral roughly by 
V / H7 ~ 0.01, i.e. the residual influence of the Hubble con- 
stant is of the order of 1% only. Since we assume that 
the universe is flat, we have fl m = 1 — fi r — In real- 
ity, the important scale is rather the distance to the last 
scattering surface (LSS), so that we should set the lower 
limit of integration to £lss ~ 1/1090. This modifies the 
result by about 2% but does not change the conclusions. 

The strength of the topological signal depends (among 
other things) on the ratio of the two scales. A universe 
that is much smaller than the causal domain leaves a 
strong imprint in the CMB, while one that is much larger 
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will be impossible to distinguish from a truly infinite uni- 
verse. The value of D corresponds roughly to the largest 
X that we can hope to detect. We show in Figures Q] [2] 
and [3] how the topology becomes easier to detect as the 
fundamental domain becomes smaller (for a fixed size of 
the particle horizon). 

For a fixed size of the fundamental domain, the pa- 
rameter affecting the ratio of the two scales is J7a as it 
changes D, see Fig. [4] If the matter density is very large, 
Q rn ~ 1, then the size of causal domain corresponds 
roughly to the size of a T [4,4,4] universe (in a purely 
matter-dominated universe, the horizon size equals twice 
the Hubble radius, hence D = 4). On the other hand, 
the causal domain of a flat low-density universe can have 
an over 1000 times larger volume. 




FIG. 4: Diameter of the particle horizon, D — 2aor/oHo/c, 
for a flat universe and as a function of J7a. This corresponds 
roughly to the maximal fundamental domain size X that we 
can hope to detect. 



The importance of the correct orientation of the fun- 
damental domain was already discussed in Ref. @. We 
found that the orientation is one of the critical parame- 
ters and needs to be taken into account in any analysis. 
This can be done by averaging quantities over the orien- 
tation (as is done for a Bayesian study), or by taking the 
maximum over all orientations (with the matched filter, 
or also when looking for correlated circles). 



on the time derivative of the gravitational potential. The 
anisotropies due to the ISW effect thus depend on the 
parameters of the background cosmology and are also 
tightly coupled to the clustering and the spatial and tem- 
poral evolution of the intervening structures. The most 
important contribution for topological studies is the late 
ISW effect due to the decay of the gravitational poten- 
tial [181 [19( . This happens in a low matter density uni- 
verse at the onset of dark energy (or spatial curvature) 
domination, when the increased rate of expansion of the 
universe reduces the amplitude of gravitational potential. 

The power spectrum of the ISW can be written as a 
function of the power spectrum of the potential P$(k) 
like 



Ct = (4tt) 2 / P^(k)k 2 dk 



2F(rj)ji(kr) drj 

' (19) 

where the LSS subscript stands for the epoch of last scat- 
tering, ji(x) is the spherical Bessel function and r is the 
comoving distance. F(k, r/) = D/x is the growth rate of 
the potential, where x is the (normalized) scale factor and 
D is the linear growth factor. The potential is related to 
the matter density contrast 6 via the Poisson equation, 
k 2 & = — |a§il m i?Q(5/x. In the matter-dominated era, 
the density contrast grows as x and therefore the gravi- 
tational potential is constant. This is no longer the case 
in the dark-energy dominated era. 

The ISW effect is seen mainly in the lowest £- values of 
the power spectrum where the cosmic variance is large. 
This makes its direct measurement very difficult. Since 
the time evolution of the potential that gives rise to the 
ISW effect may also be probed by observations of the 
large scale structure, one expects the ISW to be corre- 
lated with tracers of the large scale structure (see, e.g. 
Ref. [10|)- This could conceivably help to remove it at 
least partially. We show in Fig. [5] the information con- 
tent H of our standard T[4,4,4] model (dashed blue line) 
and compare it to a fictitious case where the ISW effect 
was turned off (red solid line; the anisotropies are then 
only generated at last scattering). The ISW is clearly 
a very strong effect that needs to be taken into account 
correctly. We also checked its influence on the detection 
algorithms discussed in section IIIBI and found the same 
result. 



Reionisation 



B. ISW effect 

The CMB photons are influenced by the change in the 
gravitational potential through which they must pass. 
The effects arising from time-variable metric perturba- 
tions are generally known as the integrated Sachs- Wolfe 
effect [TH in the linear regime and go by the names of 
Rees-Sciama effect [l6| and the effect from moving ha- 
los pjj in the non-linear regime. The ISW effect depends 



At reionisation, CMB photons scatter off free electrons, 
altering the CMB temperature anisotropies and induc- 
ing a polarised signal. The relevant effects in the case 
of topological studies are mostly those induced by reion- 
isation on the CMB temperature anisotropies. At reion- 
isation, the photon-electron interactions randomise the 
directions of propagation of the fraction of rescattered 
photons. As a result, the primary CMB anisotropies are 
suppressed and the power in the acoustic peaks is damped 



6 




10 20 



FIG. 5: Influence of the ISW effect on the information con- 
tent: The dashed blue line shows the H for the T[4,4,4] model 
while the red line shows the information content for the same 
model but without ISW effect. The difference is very impor- 
tant. 



by a factor exp (— 2r), where r is the Thomson optical 
depth. Scales approaching the horizon at reionisation 
are barely affected by the scattering, the suppression is 
essentially efficient at small scales. 

We investigate the effect of reionisation on the de- 
tectability of a non-trivial topology in Fig. [5] by com- 
paring a universe with no reionisation (red solid line) 
and a universe with an optical depth compatible with 
the WMAP1 constraints (r = 0.166, blue dashed line). 
The induced changes are negligible. The latest WMAP 
results suggest a significantly lower optical depth, there- 
fore, the effect of reionisation is even smaller. 

D. Tensor modes and defects 

Both gravitational waves and defects induce pertur- 
bations not primarily at decoupling, but throughout the 
whole evolution of the universe. Much of their contri- 
bution to the temperature anisotropies is therefore due 
to their late-time behaviour. In this they behave simi- 
larly to the integrated Sachs- Wolfe effect, and their con- 
tribution to the ai m will be mostly uncorrelated with the 
topological signal. 

In a topologically non-trivial universe defects can 
stretch across the full fundamental domain (if they are 
at least one-dimensional). This may lead to a preferen- 
tial alignment in certain directions. It also means that 
they cannot be larger than the size of the fundamental 
domain, so that their scaling behaviour has a cut-off at 
that scale. When the universe becomes larger, the defects 
will generically disappear (up to a possible net winding 
charge), leading to less fluctuation power on large scales. 
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FIG. 6: Influence of reionisation on the information content: 
The dashed blue line shows the H for the T[4,4,4] model while 
the red line shows the information content for the same model 
but without reionisation. The difference is negligible. 



Analyses [U [22j have restricted the contributions of 
primordial tensor modes and defects to be subdominant, 
and future measurements will improve limits substan- 
tially. We will consider them as a negligible source of 
"noise" for measurements of the topology. If ever their 
presence is detected, we should take them into account 
explicitely. 



IV. SENSITIVITY TO COSMOLOGICAL 
PARAMETERS 

We choose a fiducial model motivated by the 
WMAP3+SNLS results for a flat ACDM model [H], with 
H = 72.4 km s" 1 Mpc" 1 , n s = 0.95, Ov = 0.75, 
il b h 2 = 0.0234, and r = 0.085. We then focus on the 
most relevant cosmological parameters and we explore 
their effects on the detection of the topological signal 
by varying in turn Hq, n$ and Q\ (keeping all other 
parameters fixed). Table U shows the range of varia- 
tion, corresponding to roughly the two-sigma limits of 
WMAP3 (with the exception of J7a where we also inves- 
tigate other values). We do not expect the normalisation 
As or the baryon density fl^h 2 to play an important role. 
In section IIII CI we further saw that there is only a small 
difference between r = and the fiducial case. We con- 
clude that variations of r are also uncritical for measuring 
topology. 

We limit ourselves to showing the results for the 
T[4,4,4] case since it is more realistic than T[2,2,2] while 
still giving a clear signal (as opposed to T[6,6,6]). The 
basic effects due to the parameter uncertainties do not 
depend on the size of the fundamental domain. We show 
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h = {0.678,0.724,0.770} 
n s = {0.916,0.95.0,0.982} 

Q A = {0.50, 0.70, 0.75, 0.80, 0.90} 40 
TABLE I: Parameter values investigated. 



in the top panel of Fig. [7] how the information content 
H depends on the three values of ns (cf table [TJ) . The 
other two panels show how such a variation of ng affects 
the detection of a T [4,4,4] topology based on our fiducial 
model. To this end we simulated maps with the differ- 
ent values of ns (and all other parameters fixed to the 
fiducial value) and then tested them against our fiducial 
model, applying the algorithms of section III Bl We find 
that the information content does not vary significantly. 
Also the detections measured by S/N (middle panel) and 
the Bayesian evidence (bottom panel) do not change for 
f max < 10 with only small changes at higher resolutions 
(within the statistical fluctuations induced by the map 
realisations). We conclude that ns can safely be fixed to 
the central value in table Q] 

The reduced Hubble parameter h intervenes both in 
the size of the particle horizon (section IIII Ap and in the 
ISW effect (section [ill B[) . We therefore expect some de- 
pendence on its value. Repeating the same analyses as 
for ns- We find that the results depend even less on h 
than on ns , see Fig. [5] Clearly, h can also be fixed to the 
central value in table [TJ 

We follow the same procedure as above, varying fl\ 
only, for the values given in table UJ In the top panel of 
Fig. [5]we plot the information content for these models, 
with the largest £l\ leading to the highest H. The larger 
CI A, the larger the causal domain, but also the bigger 
the ISW contribution. The graph suggests that the ef- 
fect due to the size of the fundamental domain dominates 
over the increase in "ISW noise" . As before, we then in- 
vestigate how the detection changes if map and template 
have a different Q\. In contrast to the two previous pa- 
rameters, we find that there is a strong variation in the 
detectability. We only have a high S/N if map and tem- 
plate agree (black solid line in the middle panel). The 
maps with the next higher (red dashed line) and lower 
(blue dotted line) values of f^A are already below a signif- 
icant detection limit, and the other cases lie even lower. 
The evidence shows the same behaviour, with the fidu- 
cial template leading to a much higher model probability 
than the others. A full analysis requires to explore all 
acceptable values of the cosmological parameters. Given 
that the topological signal remains nearly constant for all 
parameters except J7a , it suffices to vary only £!a within 
the 2cr region allowed by WMAP3+SNLS. 

V. CONCLUSIONS 

We have studied cubic torus universes, investigating 
the change of the topological signal with respect to the 
following effects: ISW, reionisation, size of causal hori- 




max 



FIG. 7: For T[4,4,4], the effect of varying ns: The top panel 
shows the information content, for the fiducial model (solid 
black line, ns = 0.950) as well as ns = 0.916 (dashed red line) 
and n s = 0.982 (dotted blue line). The middle panel (S/N) 
and bottom panel (evidence) show how the detection changes 
when testing simulated maps against the fiducial model. The 
black solid line uses a map derived from the fiducial model, 
while the red dashed and blue dotted lines are for the smaller 
and larger ns respectively. 



zon D, topological defects and primordial gravitational 
waves. We consider these to be the potentially most rel- 
evant effects for studying topology. We find that reioni- 
sation has negligible effect. The ratio D to the size of the 
fundamental domain is the key quantity which limits our 
ability to detect a given topology. The ISW effect in turn 
is not correlated with the topological signal and so acts 
as an additional noise. The same is true for topological 
defects, and tensor modes, but they are experimentally 
constrained to a lower amplitude than the ISW effect. 
We conclude that only the ISW effect and the size of the 
causal horizon have to be considered. 

An additional important parameter is the relative ori- 
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are required for a complete analysis. 





FIG. 8: The effect of varying h for T[4,4,4]: The top panel 
shows the information content, for the fiducial model (solid 
black line, h = 0.724) as well as h — 0.678 (dashed red line) 
and h = 0.770 (dotted blue line). The middle panel (S/N) 
and bottom panel (evidence) show how the detection changes 
when testing simulated maps against the fiducial model. The 
black solid line uses a map derived from the fiducial model, 
while the red dashed and blue dotted lines are for the smaller 
and larger h respectively. 



entation of the fundamental domain Q. The methods 
we use take this already into account so that we do not 
discuss it in detail. 

When searching for traces of a non-trivial topology, 
we have to test a priori templates for all cosmological 
models, varying the full set of cosmological parameters. 
We have shown here that, given the constraints from 
WMAP3+SNLS, only the allowed variation of Q A is able 
to affect significantly the detections. This parameter al- 
ters both the magnitude of the 1SW effect and the size 
of the causal horizon. We find that the latter effect is 
more important and templates for several values of Q\ 



FIG. 9: The effect of varying fi A for T[4,4,4]: The top panel 
shows the information content, for the fiducial model (solid 
black line, Qa = 0.75) as well as Oa = 0.80 (dashed red 
line), Qa = 0.90 (dot-dashed magenta line), S1 A = 0.7 (dot- 
ted blue line), and 57a = 0.50 (long-dashed cyan line). The 
middle panel (S/N) and bottom panel (evidence) show how 
the detection changes when testing simulated maps against 
the fiducial model. The black solid line uses a map derived 
from the fiducial model. The line styles of the other models 
are the same as those of the information content (top panel) . 
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